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Topological Analysis of Computed Three-Dimensional Viscous 
Flow Fields 


GEORGE S. DEIWERT 

NASA Ames Research Center 
Moffett Field, California 


Summary 

Computed solutions of the time-dependent, Reynolds-averaged 
Navier-Stokes equations for three-dimensional flows having 
thin shear layers are analyzed, using topological concepts. 
Specific examples include the transonic flow over a body of 
revolution with conical afterbody at moderate angles of inci- 
dence to the free stream. Experimental flow-visualization 
techniques are simulated graphically to visualize the computed 
flow. Scalar and vector fluid dynamic properties, such as 
pressure, shear stress, and vorticity on the body surface, are 
presented as topological maps, and their relationship to one 
another in terms of orientation and singular points is dis- 
cussed. The extrapolation from these surface topologies 
toward the understanding of external flow-field behavior is 
discussed and demonstrated. 

Introduction 

Computational fluid dynamics has progressed in recent years 
such that with present numerical algorithms and modern compu- 
tational facilities it is becoming feasible to simulate 
complex viscous flows over three-dimensional configurations of 
practical interest. Typical solutions are obtained in the 
form of very large data sets, of the order of millions of 
words, that define the flow-field properties at a large number 
of discrete points on and about the body of interest. To 
visualize flow-field features, techniques devised and used by 
experimentalists can be emulated, and graphic displays can be 
generated that yield comparable information. Experimental 
visualization techniques include surface-oil-streak patterns, 
shadowgraphs, Schlieren photographs, interferograms , holograms, 
tufts, dye, smoke, vapor screens, and hydrogen bubbles. In 
addition, flow fields are probed for quantitative information 
with tubes, wires, films, and light beams. Similar quantitative 
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information can be extracted from computed data bases in a 
straightforward manner. 

Particular information concerning complex three-dimensional 
flow features is found in the flow-field topology on or very 
near the surface. This topology is the foundation upon which 
the external flow rests, and an understanding of the surface 
features can provide substantial insight into the behavior of 
the rest of the flow field. The analysis by Squire [l] shows 
that surface oil-flow visualizations closely approximate the 
skin-friction lines and limiting streamlines. Some considerable 
analyses by Legendre [2], Lighthill [3], Perry and Fairlie [4], 
Hunt et al. [5], and Tobak and Peake [6] provide a mathematical 
basis that permits detailed descriptions of three-dimensional 
separation and reattachment patterns based on surface vector 
patterns and their singular points. By extracting surface or 
near-surface vector and scalar topologies from discretized 
computed three-dimensional flow fields, the computational fluid 
dynamicist can use these mathematical concepts to analyze 
regions of separation and reattachment and to relate computed 
patterns to experimental observations. 

To illustrate some of the topological principles for computed 
data bases and to relate the behavior of various vector and 
scalar topological patterns to each other, a series of computed 
results for flow over a conical afterbody at moderate angles 
of incidence to the free stream is considered. 

The thin-shear-layer approximation to the time-dependent, 
three-dimensional, compressible, Reynolds-averaged Navier- 
Stokes equations was solved numerically, using the implicit 
finite-difference method of Beam and Warming [7]. The compu- 
tations correspond to a series of experiments by Shrewsbury [8] 
for flow over a pointed cylindrical forebody of caliber nine 
followed by a 15° conical afterbody and a two-thirds-forebody- 
diameter cylindrical sting. The flow conditions correspond to 
a free-stream Mach number of 0.9 and a Reynolds number based 
on forebody diameter of nearly 3 x 10 6 . Steady-state solutions 
were computed for angles of incidence to the free stream of 
0°, 2°, 4°, 6°, 8°, and 12°. Quantitative comparisons of 
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computed and experimentally measured surface pressure data for 
angles of incidence up to 6° have been reported by Deiwert [9] 
and good agreement was found. Computed results for the 8° and 
12° cases will be used as illustrative examples here. 

Discussion 

To relate various surface properties, the incompressible 
Nayier-Stokes equations are written for the surface where all 
velocity components are zero and = 0: 

P 5 = UU ?? = UtoVc 

and 


P 


n 


y v 




p( Vc • 


where E, and n are coordinates in the plane of the body 
surface and z, is normal to the surface. The above relations 
tell us that the tangential vorticity is aligned in the direc- 
tion of surface isobars. For a small distance £ from the 
surface, the velocity vector field in the plane parallel to 
the surface is given by 

V = (w w x n)c = r w c/y 

->■ . 
where co is the vorticity vector at the surface and n is 

w 

the unit surface normal vector. This relation tells us that 

-> ( -y 

(surface shear) is orthogonal to oj (tangential vorticity) 
and hence to the direction of t> e surface isobars. 

Discretized solutions obtained from finite-difference equa- 
tions do not provide detailed descriptions of behavior at 
singular points. Small-perturbation analyses of the differ- 
ential equations in the vicinity of singular points can pro- 
vide more detail and enhance our understanding of the flow 
behavior in these critical regions, which occur when the flow 
separates and reattaches. Two types of separation typically 
occur in three-dimensional flows: one when the components 

of tangential vorticity approach zero simultaneously, and the 
other when there is a convergence of surface streamlines. In 
the latter type of separation we can see from continuity that 
for flow through a rectangular stream tube whose .base b, is 
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the distance between two adjacent limiting streamlines and 
whose variable height is £ , 

1 to Is 2 b/2 = constant . 

1 W ' 

Thus, when surface streamlines converge, b decreases and 

there is a corresponding increase in X, and a "piling up" of 

stream tubes that results in three-dimensional separation in 

the form of a vortex sheet. In the former type separation 

(a i -* 0) , small-perturbation analysis suggests that in the 
w 

vicinity of a singular point there is also a pressure extremum 
which can be defined by p = p = 0, where C and n are 

~ t, n 

mutually orthogonal. 

For the present cases, limiting streamlines were determined 
by tracing particle paths parallel to the body surface (neg- 
lecting the surface normal velocity component) at a very small 
distance from the surface (c/d = 0.0000 8). These lines are 
identical to skin-friction lines and, according to Squire [l], 
should closely approximate oil-streak patterns. Shown in 
Fig. 1 are side views of limiting streamline patterns over 
the complete configuration at angles of incidence of 8° 

(Fig. la) and 12° (Fig. lb) . In both cases the streamlines 
on the cylindrical body strongly converge near the generator, 
45° from the lee generator, indicating a three-dimensional 
separation in the form of a vortex sheet. 

Immediately aft of the cylinderical-f orebody/conical-af terbody 
junction, the streamlines turn abruptly and completely circum- 
scribe the afterbody. This again is an indication of three- 
dimensional separation and will be described in more detail 
subsequently. Just aft of the conical-afterbody/cylindrical- 
sting junction, the flow reattaches along a line circumscrib- 
ing the sting, and the streamlines proceed to converge down- 
stream along a lee generator of the sting. 

Shown in Fig. 2 are surface topologies of several fluid 
dynamic properties for the afterbody at 8° incidence to the 
free stream. Included are lee, lateral, and windward views 
showing lines of surface shear (solid) with lines for p^ = 0 
and p^ = 0 (dashed) superposed, where £ and n -are the axial 
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and circumferential directions, respectively. Just aft of the 
forebody/af terbody junction there is a line of separation that 
completely circumscribes the afterbody. This line begins at a 
separation saddle point on the lee generator and terminates at 
a separation node on the windward generator. Immediately up- 
stream of this separation line is a line of p^ = 0 (dashed) 
which intersects P n = 0 lines coincident with the lee and 
windward generators and p^ = 0 lines a few degrees to either 
side of each generator. All saddle and node points are de- 
noted by symbols and are defined in the figure captions. On 
the lee side, the pressure extrema are of the form saddle-node- 
saddle and correspond to the separation saddle point indicated 
by the limiting streamlines on the lee generator. On the 
windward side, the pressure extrema are of the form node-saddle- 
node and correspond to the separation node indicated by the 
limiting streamlines on the wind generator. Aft of the 
afterbody/sting junction there is a line of reattachment com- 
pletely circumscribing the sting. There are reattachment 
nodes on both the lee and windward generators and a reattach- 
ment saddle 60° from the lee generator, as indicated by the 
limiting streamlines. Immediately downstream of the line of 
reattachment is the dashed line of p^ = 0, which intersects 
Pji — 0 lines coincident with the lee and windward generators 
and the p^ = 0 line at 83° from the lee generator. These 
points are node, node, and saddle, respectively, and corre- 
spond to the reattachment singular points directly. Between 
the line of separation and the line of reattachment, the flow 
over the surface is reversed and a separation bubble exists. 

The directions of the limiting streamlines under the bubble 
indicate counterclockwise circulation when seen in the lateral 
view. Downstream of reattachment there is a strong conver- 
gence of streamlines about 55° from the lee generator, indi- 
cating separation of the second kind. 

Shown in Fig. 3 is a lateral perspective of the afterbody at 
8° incidence. Figure 3a is a surface pressure map, in which 
C p is shown as a function of axial distance (?) and circum- 
ferential angle (cj>) from the lee generator. Surface pressure 
contours (dashed) are mapped on the afterbody in Fig. 3b and 
are superposed over the limiting streamlines to illustrate 
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regions of orthogonality. The sonic surface envelope over the 
afterbody is shown in Fig. 3c, with skin-friction lines drawn 
on the afterbody surface. It is the termination of the super- 
sonic flow inside the sonic surface by means of a three- 
dimensional shock front that induces flow separation on the 
afterbody near the forebody/afterbody junction. Included also 
in Fig. 3c is a dashed line for (to • n) =0, which is aligned 
with the vortex sheet that separates from the sting. There 
are also zero vorticity lines near both the lee and windward 
generators that cross these generators at the singular points 
indicated in the limiting streamline patterns. 

As the angle of incidence is increased, the complexity of the 
surface topology can increase. Shown in Fig. 4 are lee, 
lateral, and windward views for the afterbody at 12° incidence. 
As observed in the 8° solution, there is a line of separation 
aft of the forebody /afterbody junction and a line of reattach- 
ment aft of the afterbody/sting junction. Both lines com- 
pletely circumscribe the body. Beginning at the intersection 
of the separation line and the lee generator and progressing 
along the separation line to the windward generator, there is 
a series of singular points as follows: node-saddle-node- 

saddle-node-saddle at <|> = 0°, 22°, 67°, 102°, 160°, and 180°, 
respectively. These points are preceded in the streamwise 
direction by pressure extrema at <J> = 0°, 30°, 145°, and 180°, 
which are node-saddle-node-saddle points, respectively. It 
is of particular interest to note that there exist more singu- 
lar points on the line of separation than there are pressure 
extrema in the vicinity. The "extra" points appear to be the 
separation node and saddle at <f> = 67° and 102°, respectively. 

In addition, the flow pattern under the separation bubble is 
broken into several distinct counterrotating zones, with the 
most prominent discriminator being a line connecting the sep- 
aration saddle at <f> = 102° and the reattachment saddle at 
<j> = 75°. Along this line there is strong convergence of the 
limiting streamlines, indicating a prominent bubble bifurca- 
tion by means of a vortex sheet. The complex three-dimensional 
flow behavior in the separation bubble obviates a simpler 
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surface analysis and a one-to-one correspondence between pres- 
sure extrema and singular points in the surface pattern. 

Along the line of reattachment there are nodal points on both 
the lee and windward generators and a saddle 75° from the lee. 
These are accompanied by pressure extrema slightly downstream 
with nodal points on both the lee and windward generators and 
a saddle-node-saddle cluster at cf> = 65° , 75°, and 97° , respec- 
tively, near the reattachment saddle. Again, as in the 8° 
case, convergence of streamlines on the sting near the <f> = 70° 
generator indicates the separation of a vortex sheet. 

Shown in Fig. 5 are perspective views. of the 12° solution. 

Figure 5a shows the surface pressure map, C ( 5 ,$). Dashed 

P 

surface pressure contours are superposed on the surface 
streamline pattern in Fig. 5b. The sonic surface and surface 
streamlines are included in Fig. 5c, with a dashed line for 
(w • n) = 0 superposed. This zero vorticity line is co- 
aligned with the location of converged surface streamlines. 
Other zero vorticity lines appear near the lee and windward 
generators, crossing these lines at the separation and re- 
attachment saddles and nodes. 

Concluding Remarks 

In the examples discussed above, we have seen that in regions 
of three-dimensional separation where the surface vorticity 

-y ~y 

goes to zero (to • n. 0) , the singular points defining the 
lines of separation are preceded by accompanying pressure 
extrema of like character (i.e., nodernode and saddle : saddle) 
and singular points defining lines of reattachment precede 
accompanying pressure extrema of like character (note that 
clusters of singular points behave as a singular point so 
that a node-saddle-node behaves as a node and a saddle-node- 
saddle behaves as a saddle). For the 8° case, there is a 
direct correspondence between the pressure extrema and the 
singular points of separation and reattachment. For the 12° 
case, where the flow pattern is more complex, more singular 
points were observed along the separation line than were 
suggested by accompanying pressure extrema. In both cases 
(8° and 12°) singular points are indicated both by the 
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limiting streamlines and by zero surface vorticity. For the 
separated vortex sheet developing from the reattachment saddle 
there is an accompanying line of zero surface vorticity co- 
aligned with the converging limiting streamlines. 
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Fig. 1. Limiting streamline pattern on complete afterbody 
model at incidence (a) to the free stream (lateral view) 

M = 0.9, Re d = 3 x 10 6 : (a) a = 8°; (b) a = 12°. 
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Fig. 2. Surface shear ( — ) 
and zero pressure gradient 
( — ) topology for conical 
afterbody at 8° incidence, 
M ro = 0.9: (a) lee; (b) 

lateral; (c) windward side. 


Fig. 3. Lateral perspective 
of flow features for conical 
afterbody at 8° incidence, 

M ro = 0.9: (a) surface pres-, 

sure; (b) isobars ( — ) and 
surface shear ( — ) ; (c) sonic 
surface, surface shear ( — ) 
and zero vorticity 
(uL • n = 0) ( — ) . 
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Fig. 4. Surface shear ( — ) 
and zero pressure gradient 
topology ( — ) , for conical 
afterbody at 12° incidence. 
Moo = 0.9: (a) lee? (b) 

lateral; (c) windward side. 


Fig. 5. Lateral perspective 
of flow features for conical 
afterbody at 12° incidence, 

= 0.9: (a) surface pressure 

(b) isobars ( — ) and surface 
shear ( — ) ; (c) sonic surface, 

surface shear ( — ) and zero 
vorticity (t % • n = 0) ( — ) . 
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